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ABSTRACT 

In  ice-covered  waters  the  usefulness  of  hydrophones  is  limited  by  the 
necessity  to  create  or  to  locate  openings  through  which  the  sensors  can 
'enter  the  water.  But  ability  to  detect  submarines  in  Arctic  waters  is  im¬ 
portant.  For  this  reason  the  effectiveness  of  an  on-ice  sensor,  namely, 
the  geophone,  has  been  studied  analytically  toward  the  assessment  of  its 
feasibility.  Comparisons  were  made  between  meager  measurements  of  ice 
noise  reported  in  the  literature  and  signal  levels  calculated  from  normal 
mode  theory.  The  studies  indicate  that  a  submarine  emitting  one-watt  or 

more  at  7  or  8  Hz,  a  frequency  range  of  low  natural  noise,  might  be  detect¬ 
able  at  a  range  of  1  km  and  possibly  as  large  as  30  km. 

Noise  data  is  lacking  particularly  for  shallow  water  areas,  where  the 

geophone  seems  most  useful,  and  experimental  measurements  are  recommended 
using  three  component  geophones. 

No  matter  what  the  signal-to-noise  ratio,  array  processing  can  be  used 
to  enhance  the  probability  of  recognizing  a  submarine.  However,  the  acoustic 
environment  causes  dispersive  waves  and  this  would  add  to  the  operational 
difficulties. 
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NCaiENCLATURE 

e  voltage 
y  displacement 
(U  circular  frequency 
natural  frequency 

h  damping  factor  of  a  geophone  or  a  layer  thickness 
h^  thickness  of  i-th  layer 
G  intrinsic  sensitivity  of  a  geophone 
k  wave  number 

a  compressional  wave  velocity  in  the  i-th  layer  of 
a  multilayered  acoustical  environment 

m  mode  number,  an  integer 

n  unit  vector  in  the  direct ioft  of  propagation  for 
the  m-th  mode 

z  unit  vect.or  in  the  vertical  direction 

k  projection  of  n  k  in  the  vertical  direction 
mz  m 

' 

k  projection  of  n  k  in  the  horizontal  direction,  the 
mr  .  m 

eigen  wavenumber 

complex  valued  eigen  wavenumber 

e  imaginary  part  of  h 
mr  mr 

vertical  component  of  particle  velocity 
horizontal  component  of  particle  velocity 
R  horizontal  distance  between  source  and  receiver 
Z  depth  of  receiver 
z'  depth  of  source 
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t  time 

M  maximum  mode  that  may  exist  for  a  given  frequency 

D  radial  term  in  the  expression  for  particle  velocity, 
m-th  mode 

F  vertical  term  in  the  expression  for  particle  veloiily, 
m-th  mode 

W  excitation  function  for  the  m-th  mode 
m 

P  power  output  of  the  source 
density  of  the  i-th  layer 

L  a  function  that  relates  W  to  the  physical  layering 
of  an  inhomogeneous  acoustical  environment,  m-th  mode 

c  phase  velocity  of  the  n-th  mode 
m 

c  phase  velocity  in  general 

phase  change  occurring  for  an  up-going  wove  reflected 
at  an  interface 

phase  change  occurring  for  a  down-going  wave  reflected 
at  an  interface 

value  of  the  eigen  wave  number  in  the  i-th  layer 

A  amplitude  of  the  m-th  mode  before  reflection 
m 

k'  amplitude  of  the  m-th  mode  after  reflection 
m 

a  standard  deviation  of  a  random  surface  that  has  a 
Gaussian  distribution  function 

£(x)  height  of  the  random  surface  at  point  x 

n  number  of  reflections  in  distance  R 

9  angle  of  incidence  of  m-th  mode 
m 
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Introduction 


In  the  Arctic  and  Subarctic  regions  an  ice  cover  would  interfere  with 
the  use  of  hydrophones  for  the  detection  and  location  of  submarines.  To  use 
air-dropped  ASW  hydrophone  systems  in  ice-covered  waters,  either  a  lead  must 
be  found  or  a  hole  must  be  made  in  the  ice  through  which  the  hydrophone  can 
enter  the  water.  Either  task  is  time  consuming.  Moreover,  leads  may  not 
always  be  found  in  the  area  of  interest  and,  unless  an  area  of  thin  ice  is 
present  entry,  to  the  water  may  be  prohibitively  difficult. 

To  overcome  the  disadvantages  of  hydrophone  installation  in  ice-covered 
regions,  we  have  considei’ed  the  use  of  seismometers  as  -suitable  on-ice  sensors 
(Goldstein  1966).  The  modern  seismometer,  or  geophone  as  the  small,  high- 
frequency  types  are  called,  produces  a  voltage  output  whose  amplitude  depends 
on  the  particle  motion  of  the  medium  to  which  the  transducer  is  coupled.  Be- 
cau.se  the  acoustical  energy  from  an  underwater  sound  source  is  transmitted  as 
a  disturbance  in  particle  motion  of  the  surrounding  media,  sensitive  geophones 
located  on  the  ice  near  the  source  should  sense  the  variations  of  ice  motion. 

The  modern  geophone,  a  moving-coil  transducer,  is  particularly  sensitive 
to  motion  in  the  1  to  100  Hz  range — the  bandwidth  in  which  most  submersible 
craft  radiate  a  good  portion  of  acoustical  energy.  These  transducers  are 
rugged,  dependable  and  fairly  low  in  cost,  the  cost  increasing  with  sensitivity 
of  the  device.  Their  low  cost  allows  them  to  be  considered  as  expendable 
items.  Furthermore,  their  ruggedness  has  already  been  exploited  in  the  design 
of  an  air-droppable  seismic  surveillance  device  manufactured  for  military 
purposes  by  R-I  Controls,  Inc.  of  Minneapolis,  Minnesota. 

The  effectiveness  of  on-ice  transducers  for  ASW  work  has  never  been 
evaluated  by  experiment.  However,  a  limited  effort  by  geophysicists  has  been 
directed  toward  the  measurement  of  natural  ice  vibration  in  the  Arctic  Ocean. 
The  data  are  relatively  meager,  but  meager  as  they  are  they  provide  us  with 
basic  input  for  an  anlytical  study. 

In  the  analytical  study  reported  here  we  have  attempted  to  estimate 
quantitatively  the  signal-to-noise  ratio  resulting  from  a  submarine  beneath 
Arctic  ice.  The  signal  amplitudes  were  evaluated  numerically  by  normal  mode 
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theory  with  the  assumption  of  a  point  source  radiator.  Amplitudes  of  the 
noise,  or  natural  vibration,  were  taken  from  data  reported  in  the  literature. 
The  signal-to-noise  ratio  was  considered  for  two  acoustical  environments: 
deep  waters  such  as  the  central  Arctic  Ocean,  and  shallow  water  such  as  con¬ 
tinental  shelf  areas  which  are  deep  enough  for  a  submerged  fleet-class  sub¬ 
marine  but  may  have  a  seasonal  ice  cover. 


Basic  Theory  of  a  Geophone 

Most  seismometers,  except  for  a  few  special  devices,  basically  consist 
of  a  mass  suspended  by  springs  or  hinges  so  that  the  mass  will  remain  more 
or  less  fixed  as  the  frame  of  suspension  moves  with  the  motion  of  the  ground 
(or  ice  in  our  case).  Sensed  is  the  relative  motion  between  the  mass  and 
the  frame. 

Seismometers  are  manufactured  for  sensing  ground  motion  with  periods 
from  one  hour  to  periods  as  short  as  0.005  second.  For  the  detection  of 
nearby  natural  and  artificial  disturbances,  seismometers  with  natural  fre¬ 
quencies  of  1  to  30  Hz  are  used.  These  devices,  commonly  called  geophones, 
will  be  discussed  in  this  section. 

Geophones  are  small,  rugged,  and  simple  in  principle  (Fig.  1).  Wien 
coupled  to  a  dynamic  medium  the  transducer's  case  moves,  and  a  relative 
motion  is  set  up  between  a  magnet  fixed  to  the  case  and  an  internal  spring- 
supported  mass  on  which  a  conducting  coil  is  wound.  Neglecting  phase  changes, 
hysteresis,  eddy  currents,  etc.,  the  voltage  induced  in  the  coil  is  propor¬ 
tional  to  the  time  rate  at  which  the  coil  cuts  the  magnetic  lines  of  force. 

The  expression  for  voltage  so  generated  can  be  derived  from  the  equation  of 
motion  of  a  damped  linear  oscillator  and  Faraday's  law  of  induction: 


e  +  2hcOf-e  *  co^e  =  yG 

where 

e  =  voltage  output 
h  =  damping  term 

=  natural  frequency  of  the  mechanical  system 
y  =  case  displacement 
G  =  intrinsic  se.'s',  itivy 
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FIG.  1  COMMERCIAL  MOVING  COIL  SEISMOMETER 


The  dots  mean  time  derivatives.  G,  the  intrinsic  sensitivity,  is  propor¬ 
tional  to  the  number  of  turns  on  the  coil,  the  dipole  moment  of  the  magnet, 
and  the  mass.  Because  the  power  output  is  proportional  to  mass,  the  voltage 
output,  hence  sensitivity,  is  proportional  to  the  square  root  of  the  mass. 
Geophones  with  extremely  high  sensitivity  have  masses  of  a  kg  or  more. 
Damping,  h,  between  the  mass  and  the  frame  is  required  so  that  the  instrument 
does  not  oscillate  freely  when  excited. 

Two  limiting  cases  of  Eq.  1  are  of  ititerest  for  they  point  out  important 
characteristics  of  a  geophone.  When  the  natural  frequency  is  much  less  than 
the  excitation  frequency,  Eq.  1  may  be  reduced  to 
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e  =  JG  (co^  «  w) 


From  this  we  see  that  the  voltage  output  is  proportional  to  the  velocity  of 
the  case  and  thus  the  particle  velocity  of  the  medium.  In  this  situation  the 
geophone  behaves  as  a  velocity  gage.  On  the  other  hand  a  different  prime 
characteristic  is  revealed  when  the  natural  frequency  is  much  larger  than 
the  excitation  frequency.  Then  Eq-  1  may  be  approximated  by 


e  =  yG/w^  »  w) 


and  the  voltage  output,  although  small,  is  proportional  to  the  third  deriva¬ 
tive  of  particle  displacement.  Thus  for  excitation  near  frequencies  of  maxi¬ 
mum  geophone  sensitivity  the  voltage  generated  will  have  a  complex 

relation  to  ground  motion:  the  voltage  output  will  be  a  superposition  of 
particle  velocity  and  acceleration  effects. 

This  behavior  is  not  detrimental  to  most  geophone  applications.  If  we 
are  interested  in  knowing  only  whether  or  not  a  disturbance  occurs,  we  need 
not  translate  the  observed  voltage  into  absolute  quantities  of  ground  motion. 
Our  only  concern  is  to  choose  a  transducer  with  a  natural  frequency  near  the 
frequency  of  maximum  excitation  expected.  However,  calibration  tests  can  be 
made  so  that  a  conversion  between  voltage  and  velocity  (or  displacement)  is 
possible. 

Natural  Particle  Velocities  of  Arctic  Ice 

Recent  experiments  have  shown  that  Arctic  ice  is  in  continuous  horizon¬ 
tal  and  vertical  motion.  In  this  report  we  shall  consider  only  the  vertical 
particle  velocity  of  the  ice  because  it  is  the  only  component  of  velocity 
that  has  been  measured  over  long  periods  of  time.  Fortunately,  particle 
velocities  induced  in  thin  ice  from  an  underwater  low-frequency  sound  source 
are  theoretically  in  the  vertical  direction  (Goldstein,  1966). 

Figure  2  shows  the  amplitude  spectra  of  vertical  particle  velocity  of 
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FIG.  2  AMBIENT  VERTICAL  PARTICLE  VELOCITY  OF  ARCTIC  ICE 


Arctic  ice  as  determined  by  Hunkins  (1962)  and  Prentiss,  et  £l  (1965).  These 
spectra  are  plotted  in  the  conventional  units  used  by  seismologists:  millimicrons 
per  second  (mp,/sec)  where  1  mp,  =;  10“®  meters.  The  general  characteristics  of  the 
observed  spectra  are  a  decrease  of  particle  velocity  with  increasing  frequency, 
a  minimum  velocity  in  the  frequency  range  of  one  to  five  Hz,  and  a  wide  range 
between  minimum  and  maximum  noise  intensities  over  the  10"^  to  10^  Hz  bandwidth. 
The  spectra  shown  were  obtained  from  observations  made  on  ice  over  the  deeper 
portions  of  the  Arctic  Ocean  where  the  water  depth  is  2.5  to  2.8  km.  This  envi¬ 
ronment  qualification  of  the  spectra  shown  is  judged  to  be  important  because  the 
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noise  spectra  depend  on  the  physical  properties  of  the  coupled  ice-to-sea- 
water-to-sea-bottom  acoustical  system.  For  example,  Milne  and  Clark  (196d) 
reported  that  closely  spaced  resonances  in  seismic  background  noise  in  the 
band  0.1  to  2.0  Hz  were  observed  at  the  bottom  of  ice-covered  estuaries  of 
the  Canadian  Arctic  where  the  measured  water  depth  was  0.48  km.  These  reso¬ 
nances  could  be  explained  as  the  result  of  multiple  internal  reflections  of 
plane  waves  within  a  layered  medium  with  acoustical  properties  found  from 
seismic  experiments. 

During  periods  of  minimum  noise  the  particle  velocity  of  the  ice  is  as 
small  as  observed  at  a  quiet  continental  site.  This  was  first  pointed  out  by 
Prentiss,  et  £l  (1965)  who  compared  their  observations  of  Arctic  ice  displace¬ 
ments  to  the  continental  crust  displacements  given  by  Brune  and  Oliver  (1959). 
During  periods  of  maximum  disturbance  the  particle  velocity  of  Arctic  ice  will 
be  several  hundred  mp,/sec,  sufficient  to  completely  overwhelm  the  signal  from 
a  nearby  underwater  source  radiating  a  few  watts  of  acoustic  power. 

The  exact  causality  between  physical  processes  and  natural  ice  vibrations 
is  not  clear,  but  Arctic  researchers  have  been  able  to  identify  a  number  of 
general  causes  for  the  vibrations.  Foremost  is  wind  action.  Tlie  wind  is 
responsible  for  much  of  the  low  frequency  (less  than  lO"®  hz)  bobbing  of  Llie 
sea  ice  (Prentiss,  e_t  £l ,  1965),  the  noise  above  1  Hz  coming  from  impacts 
between  adjacent  floes,  and  for  the  motion  of  granular  snow  on  the  ice  surface 
(Mi.lne  and  Canton,  1964). 

Temperature  variations  also  seem  to  affect  the  level  of  natural  vibration 
amplitudes.  During  the  winter  months  the  falling  air  temperature  late  in  the 
afternoon  is  accompanied  by  an  increasing  vibration  noise  level.  Tliis  is 
attributed  to  ice  cracking  as  a  result  of  increased  tensile  stresses  (Milne 
and  Canton,  1964). 

The  particle  velocities  of  the  natural  vibration  have  a  seasonal  variation, 
the  average  velocity  amplitude  being  larger  in  the  winter  than  in  the  late 
sumner  (Milne  and  Canton,  1964;  Mellen  and  Marsh,  1965).  In  the  late  summer 
the  ice  cover  becomes  broken  and  "fluidized”,  conditions  which  make  the  ice 
a  poor  generator  of  noise  (Mellen  and  Marsh,  1965)  and  a  poor  transmission 
medium  for  ice  vibrations. 
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The  Propagation  of  Sound  in  the  Arctic  Ocean  and  Shallow  Arctic  Waters 

In  the  ice-covered  Arctic  Ocean  the  propaRation  velocity  of  acoustic 
waves  in  the  water  increases  bilinearly  with  depth  beneath  the  ice  cover. 
Figure  3  shows  the  measured  variation  in  propagation  velocity  and  density 
in  the  Arctic  Ocean.  For  the  crustal  portion,  the  data  are  based  on  seismi 
refraction  surveys  in  the  Bering  Sea. 


DENSITY  —  g /cn? 

'  TB-‘)89522-a2R 


FIG.  3  VELOCITY  AND  DENSITY  VARIATION  IN  THE  CENTRAL  ARCTIC  OCEAN 
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In  the  shallow-water  regions  of  the  Arctic  the  acoustic  velocity  of 
the  ice-covered  waters  probably  increases  slightly  with  depth  as  in  the  case 
of  the  deep  waters.  However,  along  the  continental  shelf  the  bottom  itself 
causes  a  pronounced  velocity  increase  at  shallow  depth. 

In  both  environments  there  is  a  low-velocity  channel  near  the  water 
surface  which  results  in  guided  wave  propagation  of  acoustical  disturbances. 

A  portion  of  the  energy  emitted  by  a  sound  source  is  ti-apped  in  the  low- 
velocity  channel  by  a  process  of  internal  reflections.  A  downgoing  distur¬ 
bance  in  particle  velocity  will  be  reflected  upward  at  some  depth.  At 
reflection,  the  wave  will  change  phase,  and  if  reflection  is  incomplete,  a 
loss  in  amplitude  will  also  occur.  A  similar  occurrence  will  be  experienced 
by  an  upgoing  disturbance  impinging  on  the  ice.  If  the  ice  thickness  is 
small  compared  to  the  wavelengths  of  sound  in  ice,  then  the  ice  has  little 
effect  on  the  reflection  coefficient  of  the  incident  wave.  Tliat  is,  we  may 
treat  the  water  surface  as  a  "free"  surface  and  thereby  assume  complete 
reflection  and  a  180-degree  phase  change  of  the  incident  wave.  This  assump¬ 
tion  was  shown  to  be  valid  for  ice  less  than  3  meters  thick  reflecting  fre¬ 
quencies  below  100  Hz  (Goldstein,  1966). 

To  be  precise,  two  waveguides  exist:  the  near-surface  water  channel 
and  the  ice  guide.  At  frequencies  less  than  100  Hz  the  phase  velocities  of 
free  (trapped)  vibrations  of  the  ice  are  either  smaller  (flexual  vibrations) 
or  greater  (longitudinal  vibrations)  than  the  phase  velocities  of  the  guided 
waves  in  the  water  (Press  and  Ewing,  1951).  Therefore,  there  is  no  coupling 
of  energy  between  the  water  and  the  ice  layer,  and  our  previous  remarks  that 
the  ice  can  be  ignored  in  the  propagation  of  underwater  sound  still  holds. 
However,  for  frequencies  above  100  Hz  we  eventually  must  account  for  energy 
refracted  into  the  ice.  Depending  on  the  ice  thickness  and  angle  of  incidence, 
the  refracted  energy  will  be  reradiated  in  varying  degrees  into  the  water  or 
trapped  within  the  ice  layer.  However,  because  of  absorption  and  the  irreg- 
ulai-  ice  surfaces  the  ice  layer  waveguide  is  extremely  lossy  for  high  fi-equency 
propagation.  We  can,  therefore,  ignore  the  ice  waveguide  at  source  ranges  of 
one  kilometer  and  more. 

In  this  report  we  shall  address  ourselves  to  the  propagation  of  acoustic 
waves  at  or  below  100  Hz.  Although  we  recognize  that  the  power  radiated  by 

8 


CONFIDENTIAL 


CONFIDENTIAL 


a  submarine  extends  to  frequencies  well  above  100  Hz  we  impose  this  frequency 
constraint  because  natural  ice-vibration  (noise)  measurements  are  reported 
only  at  frequencies  below  100  Hz,  and  because  most  of  the  submarine  radiated 
power  is  also  below  100  Hz.  With  this  constraint,  then,  we  outline  below 
the  basic  theory  and  equations  lor  guided  wave  propagation  as  they  apply  to 
deep  (2.8  km)  and  shallow  (0.1  km)  ice-covered  waters.  In  later  sections  we 
will  consider,  semi-quantitatively,  the  use  of  on-ice  sensors  at  frequencies 
above  100  Hz. 

For  waves  of  constant  frequency  there  will  be  only  a  finite  number 
of  ray  directions  that  lead  to  the  propagation  of  a  large  amplitude  dis¬ 
turbance.  The  necessary  condition  for  such  a  path  is  that  constructive  phase 
interference  occurs  at  all  points  between  the  upward  and  downward  reflecting 
interfaces.  Clay  (1964)  describes  this  condition  in  terms  of  a  phase  integral 

j  k(n  ‘  'z)dz  -  mV  (2) 

0 

where 

h  =  tlie  vertical  distance  between  reflecting  surfaces 
k  -  the  wave  number  =  co/d 

a  =  the  compressional  wave  velocity  in  the  layer,  assumed 
here  to  be  constant 

=  a  unit  vector  in  the  direction  of  wave  propagation  tliat 
leads  to  constructive  interference 

■f  =  a  unit  vector  in  the  vertical  direction 

m  =  an  integer,  the  mode  number 

a  nd 

k(Ti  •  t)  s  fe 

A  A  Z 

Equation  2  carries  no  distinction  between  wave  directions  that  lead  to 
complete  internal  reflections  and  those  that  do  not.  We  assume  complete 
reflection  at  the  surface  for  all  angles  of  incidence;  but  at  the  lower 
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reflecting  interface  there  will  be  a  critical  angle  of  incidence  below 
which  a  portion  of  the  energy  will  be  refracted  away  with  each  reflection. 

At  long  ranges  from  the  source  the  contribution  from  energy  propagating  at 
angles  less  than  critical  will  be  small.  Therefore,  we  will  include  in  our 
analysis  only  the  contribution  of  those  waves  unattenuated  by  bottom  refrac¬ 
tions.  With  this  constraint  there  exists  for  each  mode,  m,  a  low  frequency 
cut-off  corresponding  to  critical  angle  reflection  (Officer,  1958). 

Acoustic  waves  satisfying  Eq.  2  will  be  dispersive,  i.e.,  different 
frequencies  will  travel  with  different  phase  and  group  velocities  and, 
depending  on  the  number  of  solutions  (modes)  permissible,  a  particular 
frequency  can  have  several  different  phase  and  group  velocities.  This  dis¬ 
persion  is  caused  by  the  variation  of  velocity  with  water  depth  and  is  re¬ 
ferred  to  as  geometrical  dispersion  (Officer,  1958). 

The  phase  integral  condition  is  the  basis  for  the  theory  of  normal 
mode  propagation.  This  approach  to  the  evaluation  of  disturbance  intensities 
propagating  in  a  waveguide  is  particulai'ly  useful  because  of  the  analytical 
solutions  developed  by  Pekeris  (1948)  and  Tolstoy  (1955,  1958),  among  others. 
With  these  solutions  we  have  been  able  to  evaluate  the  surface  particle 
velocity  amplitudes  from  a  point  source  within  models  of  the  Arctic  Ocean  and 
shallow-water  Arctic  areas.  These  amplitudes  have  been  calculated  at  ranges 
for  which  the  normal  mode  approach  is  applicable,  at  horizontal  ranges  greater 
than  several  layer  thickness.  In  this  section  a  brief  exposition  of  the  ana¬ 
lytical  approach  will  be  given.  Additional  material  concerning  the  analysis 
is  found  in  Appendix  A. 

At  distances  from  the  source  for  which  unattenuated  normal  mode  propagation 
is  applicable,  the  particle  velocity  signal  at  frequency  w  will  be  a  super¬ 
position  of  signals  from  .M  modes.  Assuming  a  simple  harmonic  point  source  the 
vertical  and  radial  components  of  particle  velocity  con  be  expressed  as: 


n-il 


cot  ik^pD^(k^^R)Fjk^^Z,k^Z’)  exp 


(3a) 


L 
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\'\(R,Z,t) 


-  to  (fc  /OF  (fe  ,Z,A-„.Z'H(4fi)"’=  ^  exp  ifth  (31.) 

n«r  n  n  z  n  *  ' 

«  “  I 


who  re 


and 


V  = 

r 

P,  = 

z  = 

2'  - 


the  vertical  coinpoiieiit  of  particle  velocity 
tl>e  horizontal  component  of  particle  velocity 
the  horizontal  range  I  roin  source  to  receiver 
tl>e  depth  of  receiver 
the  depth  of  source 
tlie  eigen  wave  number 


4 


k 


2 

«  r 


The  value  M  will  be  finite  but  will  increase  with  increasing  tu*  As  we 

see  from  Eq.  3a  and  3b  the  expressions  for  particle  velocity  are  composed 

of  two  terms  D  and  F  that  have  only  radial  and  vertical  arguments,  respec- 
m  m 

tiveiy.  At  large  values  of  R  the  radial  function  can  be  approximated  by: 
D^(k^^R)  =  /?"^  exp  +  V'l)] 


which  shows  that  these  waves  spread  cylindrically  (as  R  instead  of  R“^) 
and  that  the  modes  are  displaced  in  phase  from  each  other  by  an  amount  that 
depends  on  the  range  and  the  waveguide  eigen  values,  environment 

for  which  guided  waves  may  exist  the  cylindrical  spreading  of  energy  permits 
the  long  range  propagation  of  sound.  The  phase  displacement  between  modes 
is  a  consequence  of  the  dispersive  nature  of  guided  wave  propagation. 

Near  the  surface,  or  within  the  low-velocity  near-surface  channel,  the 

function  P  can  be  expressed  as: 
m 
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sii,  (k  Z)  s 


I II 


<,.i ) 


The  sinusoidal  deptli  dependencies  of  this  function  are  duo  to  the  standing 

wave  nature  of  the  guidea  modes.  F  is  a  maximum  whenever  source  and  receiver 

m 

are  at  antinodes  and  a  minimum  at  nodes  of  the  standing  wave  pattern  dictated 

by  the  values  of  k  that  satisfy  Eq.  2.  The  waves,  of  course,  are  not  standing 
mz 

waves  as  may  be  excited  in  an  organ  pipe;  they  are  traveling  waves.  From  Eqs. 

3a,  3b,  and  5,  it  can  be  seen  that  the  vertical  particle  velocity  will  be  a 
maximum  and  the  radial  particle  velocity  will  be  zero  at  the  surface  for  a  source 
at  any  depth  in  the  near-surface  low-velocity  channel.  This  is  a  consequence  of 
the  free-surface  boundary  condition  which,  for  a  thin  ice  cover,  is  a  valid 
assumption  at  frequencies  below  100  Hz.  Because  is  zero  at  the  surface  we 
shall  ignore  the  radial  component  of  velocity  in  subsequent  discussions. 


The  function  W  (k  h)  in  Eq.  5  is  the  excitation  function  for  the  m-th 
m  mz 

mode.  It  is  a  measure  of  the  power  radiated  into  the  mode  and  depends  on  the 

source  power  output,  the  frequency  of  excitation,  and  the  physical  parameters 

of  the  medium.  A  simplified  approach  to  the  evaluation  of  the  eigen  wave  numbers 

and  Wm  (k  h)  is  first  to  anproximate  the  inhomogeneous  acoustical  environment 
mz 

by  N  horizontal  layers,  each  with  a  constant  velocity  and  density.  By  moans  of 

this  approximation  Tolstoy  (1955)  describes  how  W  may  be  numerically  evaluated 

m 

for  an  arbitrary  system  of  N  layers  (Appendix  A).  Pekeris  (1948)  gives  an 

exact  algebraic  expression  for  W  for  the  two-  and  three-layer  cases,  and 

m 

Pekeris'  solution  was  used  in  the  calculations  presented  here.  In  the  two- 

layer  case,  only  one  lov/er  reflecting  interface  at  a  depth  h  be/.eath  the  surface, 

the  function  W  is: 

m 


W  {k  b)  =  — 


2  / 


k  P  t 

nr'  1 


k  L"* 

tt  z  m 


(6) 


where 


=  k  h  -  sin  (fe  h)  cos  {k  h)  -  (p,  sin  (k  li)/pA  tan  (.k  b) 

n  z  m  z  n  z  •  c  n  z 
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Tne  a)~'  factor  occurs  because  we  have  purposely  chosen  to  represent 
the  source  by  one  whose  power  output  is  constant  at  all  frequencies  in  an 
infinite,  homogeneous  medium*.  This  was  done  so  that  the  computed  particle 
velocity  intensities  at  various  frequencies  could  be  directly  compared  for 
unit  power  output. 

Because  our  objective  is  to  obtain  a  measure  of  V  at  the  surface  (Z- ^’) 

z  it 

caused  by  the  superposition  of  all  possible  modes,  we  ignore  the  e  tei'in 
in  Eq.  3a  and  solve  for  the  modulus  of  this  complex  expression: 


il^(/?,0)|  =  (V\  • 


where  V  is  the  complex  conjugate  of  V  .  This  multiplication  yields  the 

Z  2 

following  representation  for  the  particle  velocity; 


liMrt.ol  = 


2 


I  sin 


X  ik  k  Z,  L  )■* 

nr  nr  n  n 


(7) 


Equation  7  shows  that  the  amplitude  Iv  (R,0)|  is  the  sum  of  two  terms 

"V  ^ 

both  of  which  attenuate  as  R  a  constant  term  and  a  fluctuating  term 
caused  by  the  phase  interference  of  the  various  modes.  This  equation  foi-ms 
the  basis  for  the  data  to  be  presented.  However,  before  discusing  the 
numerical  results  it  is  important  to  consider  some  of  the  restrictions  on 
the  accuracy  of  the  numerical  evaluations. 

In  general,  the  amplitude  fall-off  that  might  be  observed  experimentally 
-V 

is  greater  than  R  for  several  reasons.  For  a  shallow  water  environment, 

*A  velocity  potential  source  with  the  form  (A/uj)exp  iout  will  yield  the  de¬ 
sired  effect  (Tolstoy,  1958). 
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attenuation  of  the  guided  waves  can  be  caused  by  the  sloping  of  the  bottom, 
incomplete  reflection  due  to  shear  waves  refracted  into  the  bottom,  roughness 
of  the  bottom  and  under-ice  surfaces,  and  the  porous  nature  of  bottom  sedi¬ 
ments.  We  can  safely  ignore  thermoelastic  effects  and  volume  scattering  for 
low-frequency  waves.  For  a  deep  water  environment,  potential  causes  for 
attenuation  are  lateral  variations  in  acoustic  impedance  and  roughness  of 
the  ice  bottom. 

To  account  for  the  effects  of  attenuation  other  than  spreading  losses, 

investigators  have  taken  several  approaches.  For  example,  to  better  describe 

wave  propagation  in  shallow  water  where  bottom  losses  occur  the  eigen  wave- 

numbers  can  be  made  complex,  h  =k  +  i  e  .  Provided  the  values  e  are 

mr  mr  mr  mr 

small,  and  experiment  seems  to  confirm  this  (Ewing,  Pekeris,  and  Worzel,  1948; 

Tolstoy,  1958),  the  eigen  values,  k^^,  are  not  perturbed  by  the  addition  of 

the  imaginary  component  and  Eq.  3  is  modified  only  to  the  extent  that  takes 

on  an  exp(-e  R)  factor, 
mr 

Brekhovsikikh  (1960)  analyzed  the  condition  >.  f  imperfect  reflections  at 
one  of  the  boundaries  and  showed  that  for  large  ranges  (large  compared  to  the 
waveguide  thickness)  the  average  decay  of  signal  intensity  could  be  represented 

-3/ 

by  a  R  fall-off.  This  fall-off  is  intermediate  between  the  cylindrical 
-1/ 

law  (R  perfect  reflection)  and  the  spherical  law  (R”^,  no  reflection). 

Clay  (1964)  extended  Eckart’s  (1953)  study  of  the  reflection  of  plane 
waves  from  slightly  irregular  surfaces  and  was  able  to  show  the  effect  of  an 
irregular  reflecting  interface  on  guided  wave  propagation.  When  such  a  surface 
exists  a  portion  of  the  guided  wave  energy  may  be  scattered  in  a  direction 
that  does  not  satisfy  the  phase  integral  (Eq.  2).  The  rate  of  attenuation 
increases  with  height  of  the  surface  irregularities,  acoustic  wave  frequency, 
and  mode  number. 

In  the  following  presentations  of  calculated  signal  (particle  velocity) 
amplitudes  we  include  the  ice  roughness  attenuation  on  the  guided  waves.  The 
development  for  the  correction  term  that  must  be  incorporated  into  Eq.  7  is 
given  in  Appendix  B.  Here  we  will  explain  only  that  for  tho  attenuation  we 
assume  the  ice  bottom  to  be  a  random  surface  with  a  Gaussian  distribution  of 
amplitudes.  The  standard  correction  factor  employed  throughout  is  based  on 
a  3-meter  standard  deviation  in  the  surface  irregularities,  a  value  proposed 
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by  Mellen  and  Marsh  (1965)  fz-om  their  experiments  in  measuring  the  scattering 
of  sound  by  Arctic  sea  ice. 


Comparison  of  Signal  and  Noise  Amplitudes  on 
the  Surface  of  an  Idealized  Arctic  Ocean 

The  Arctic  Ocean  acoustical  environment  of  Fig.  3  was  partitioned  into 
three  homogeneous  layers  with  properties  shown  in  Table  1. 


Table  1 

The  Approximation  for  the  Deep  Water  Environment 


Layer 

Velocity 

(m/sec) 

Density 

(gm/cm^) 

Thickness 

(m) 

1 

1441 

1.028 

350 

2 

1480 

1.039 

2450 

3 

1525 

1.400 

OO 

The  surface  layer  was  given  a  velocity  and  density  corresponding  to  the 
average  values  found  in  the  region  of  the  near-surface  segment  of  the 
roughly  bilinear,  oceanic  velocity  variation  (Fig.  3).  To  the  middle  layer 
we  assigned  a  velocity  and  density  equal  to  the  average  values  of  the  lower 
segment  of  the  bilinear  velocity  variation.  The  bottom  layer,  or  half¬ 
space,  was  given  the  approximate  acoustical  properties  of  water-saturated 
oceanic  sediments.  It  turns  out  that  the  magnitude  of  velocity  and  density 
assigned  the  ocean-bottom  rocks  has  a  negligible  influence  on  the  normal 
mode  propagation  because  this  region  does  not  reflect  the  guided  waves.  That 
is,  as  we  have  formulated  the  problem,  the  lower  reflecting  surface  of  the 
waves  is  the  interface  between  layers  1  and  2.  This  is  a  simplification  of 
the  true  conditions  but  is  judged  adequate  for  the  objective  of  this  study. 

Shown  in  Fig.  4  are  the  phase  velocity  dispersions  for  the  first  11 

modes  calculated  for  the  simplified  acoustical  representation  of  the  Az-ctic 

Ocean.  Tlie  eigen  wave  numbers  k  are  related  to  the  phase  velocities  bv 

mr  •' 

the  I'elation  intermediate  step  toward  the  numerical  evalu¬ 

ation  of  the  particle  velocity  intensities  at  the  surface  of  the  water,  the 
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FIG.  4  CALCULATED  PHASE  VELOCITY  DISPERSION  OF  GUIDED  WAVES  IN  DEEP  ARCTIC  WATER 


excitation,  a  measure  of  the  signal  strength  of  each  mode,  was  calculated 
using  the  eigen  wave  number  information  (see,  for  example,  Eq.  3a,  6). 

These  excitation  strengths,  shown  in  Fig.  5,  represent  a  mode-by-mode  nu¬ 
merical  evaluation  of  Eq.  3a  for  a  source  emitting  a  constant  one-watt 
acoustical  power  at  all  frequencies.  In  this  representation  we  ignore  the 
time  dependency  and  set  equal  to  unity  to  eliminate  the  effect  of  range  and 
and  phase.  The  choice  of  one-watt  in  the  calculations  is  made  with  the 
realization  that  one  watt  is  probably  more  power  than  found  in  individual 
lines  of  a  submarine  line  spectrum.  However,  unit  power  output  is  used 
here  and  throughout  all  subsequent  discussions,  and  the  numerical  results 
presented  may  be  adjusted  to  any  desired  output  power  merely  by  multiplying 
by  the  square  root  of  that  power  in  watts  (Eq.  6). 

The  excitation  strengths  illustrate  the  transparency  of  the  Arctic  Ocean 
environment  for  guided  wave  propagation  with  respect  to  mode  number,  frequency, 

16 


CONFIDENTIAL 


excitation  strength 


CONFIDENTIAL 


and  source  depth..  For  example,  Fig.  5  shows  how  changing  source  depth 
affects  the  excitation  of  10  Hz  normal-mode  vibrations.  As  the  source 
depth  increases  from  25  to  125  meters  the  resulting  signal  amplitude  will 
increase  by  a  factor  of  almost  five. 
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FIG.  5  CALCULATED  RELATIVE  EXCITATION  STRENGTHS  OF  A  FEW  MODES 
FOR  VARIOUS  SOURCE  DEPTHS  IN  DEEP  ARCTIC  WATER 

Of  more  specific  interest  than  the  curves  in  Fig.  5  are  the  range- 
dependent  amplitudes  [V^l  for  particular  frequencies  calculated  according 
to  Eq.  7.  Shown  in  Fig.  6  are  the  results  of  these  calculations  for  a  25-Hz 
source  at  depths  of  25,  50,  75,  and  125  meters.  At  this  low  frequency  only 
three  modes  contribute  to  the  total  signal  amplitude,  and  comparative  calcu¬ 
lations  showed  that  ice-bottom  roughness  had  negligible  effect  on  signal 
amplitudes  out  to  the  maximum  range  considered.  Included  in  Fig.  6  for 
comparison  are  the  signal  amplitudes  calculated  for  the  case  where  the  top 
layer  becomes  infinitely  thick:  Curves  A,  B,  and  C. 

To  illustrate  the  relative  magnitude  of  the  calculated  signal  amplitudes 
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and  the  observed  ambient  noise  amplitudes,  the  maximum  and  minimum  noise 
levels  at  25  Hz  are  indicated  in  Fig.  6.  Tliese  data  were  taken  directly 
from  the  curves  in  Fig.  2.  Also  included  in  Fig.  6  is  a  decibel  scale  which 
gives  the  signal  amplitude  in  terms  of  dB's  above  the  minimum  noise  level 
observed.  A  general  conclusion  is  that  because  the  signal  amplitude  from 
the  assumed  1-watt  source  lies  close  to  the  minimum  natural  noise  level 
there  will  be  periods  when  the  background  noise  will  completely  overwhelm 
the  signal. 
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FIG.  6  VERTICAL  PARTICLE  VELOCITY  AS  A  FUNCTION  OF  RANGE 
FOR  VARIOUS  SOURCE  DEPTHS  (Three  Modes,  25  Hz) 


To  develop  a  less  general  conclusion  we  would  need  to  know  the  most 
likely  noise  amplitude  occurring  at  a  given  time  of  day  and  season  of  year. 
Existing  data  is  insufficient  to  permit  a  complete  analysis,  but  Prentiss' 
data  have  been  used  to  plot  observed  noise  amplitudes  versus  local  time, 
observed  from  26  April  through  14  May  1962  (Fig.  7).  The  data  suggest  no 
correlation  between  noise  intensity  and  time  of  day.  However,  the  data  at 
any  instant  weie  well  above  the  minimum  level.  As  an  arbitrary  dividing  point. 
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FIG.  7  AVERAGED  VERTICAL  .^ARTICLE  VELOCITIES  OF  25-Hz  NATURAL  ICE 
VIBRATION,  26  APRIL  TO  14  MAY,  (After  Prentiss  et  al,  1965) 


100  m(i/sec  was  chosen  as  the  separation  between  noisy  and  quiet  conditions.  1 

Considering  just  the  amplitudes  during  the  "quiet"  intervals  we  find  that  r 

average  noise  is  approximately  65  m(j,/sec;  a  noise  amplitude  that  is  uniformaly  |, 

greater  than  the  calculated  amplitudes  signal  except  at  very  short  ranges  from  (, 

I''; 

the  source.  Thus,  even  under  essentially  quiet  conditions  a  source  output  of  (: 


at  least  4  watts  would  be  required  to  detect  a  source  within  20  km. 

The  unfavorable  signal-to-noise  relation  worsens  when  we  consider  fre¬ 
quencies  of  50  and  100  Hz.  For  50-Hz  propagation  the  signal  amplitudes  corre¬ 
sponding  to  several  source  depths  are  plotted  against  range  in  Fig.  S.  In 
Fig.  9  it  may  be  seen  how  the  ice-roughness  attenuation  alters  the  numerical 
value  for  1 V^  |  :  scattering  of  50-Hz  acoustic  waves  by  the  ice  bottom  leads 
to  slightly  smaller  signal  amplitudes. 

For  100-Hz  propagation,  the  effect  of  ice  roughness  is  more  inhibiting 
to  the  signal  amplitudes.  Figures  10  and  11  show  that  the  ice  bottom  rough- 
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Beyond  a  range  of  20  km  the  attenuation  amounts  to  6  dB  or  more  In  places. 
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FIG.  8  VERTICAL  PARTICLE  VELOCITY  AS  A  FUNCTION  OF  RANGE 
FOR  VARIOUS  SOURCE  DEPTHS  (Six  Modes,  50  Hz,  Ice  Bottom 
Irregularities  Ignored) 
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FIG.  9  VERTICAL  PARTICLE  VELOCITY  AS  A  FUNCTION  OF  RANGE: 
SO-Hz  OUTPUT,  SOURCE  DEPTH  25  METERS  (6  Modes) 
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FIG.  10  VERTICAL  PARTICLE  VELOCITY  AS  A  FUNCTION  OF  RANGE: 
lOO-Hz  OUTPUT,  SOURCE  DEPTH  25  METERS  (11  Modes) 
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FIG.  11  VERTICAL  PARTICLE  VELOCITY  AS  A  FUNCTION  OF  RANGE: 
100-Hz  OUTPUT,  SOURCE  DEPTH  50  METERS  (11  Modes) 
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At  both  50  and  100  Hz  the  signal  intensities  remain,  in  general,  above 
the  observed  minimum  noise  level  out  to  the  maximum  ranges  considered,  but 
the  signal  intensities  are  so  close  to  the  minimum  noise  level  that  unless  a 
condition  of  extreme  quiet  exists  the  weak  signal  would  not  be  perceptible 
above  the  ambient  noise.  Again  using  the  information  given  by  Prentiss  e_t  al^ 
(1965)  we  find  the  average  noise  was  235  mjx/sec  at  50  Hz  for  six  recording 
sessions  between  18  and  27  April  1962.  During  the  same  time  interval  the 
average  noise  was  257  mjji/sec  at  100  Hz,  based  on  only  four  recording  sessions. 
Although  the  data  are  too  few  to  yield  a  reliable  estimate,  the  average  noise 
given  suggest  that  unless  the  source  emits  far  more  power  than  1  watt  at  100  Hz 
its  detection  will  be  hampered  by  the  presence  of  ambient  noise. 

The  best  signal-to-noise  ratio  exists  at  relatively  low  frequencies,  5  to 
10  Hz.  In  Fig.  2  the  maximum  observed  noise  level  is  shown  to  drop  sharply  at 
10  Hz,  falling  from  about  200  m)j,/sec  at  10  Hz  to  10  mp,/sec  at  7  Hz.  As  guided 
wave  propagation  is  strongly  excited  at  7  Hz  (Fig.  5)  it  seems  that  an  on-ice 
transducer  should  be  most  useful  if  the  source  sought  emits  strong  spectral 
lines  at  7  or  8  Hz.  A  source  emitting  a  one-watt  line  at  7  Hz  will  yield  a 
particle  velocity  signal  at  short  ranges  (<  30  km)  that  we  estimate  from  the 
curves  to  be  roughly  equal  to  the  maximum  noise  level. 

Comparison  of  Signal  and  Noise  Intensities  on  the 
Surface  of  an  Idealized  Shallow  Water  Arctic  Region 

There  exists  in  the  A.rctic  and  Subarctic  a  large  total  area  of  relatively 
shallov/  water  which  has  a  seasonal  cover  of  land-fast  or  pack  ice.  Because 
parts  of  these  shallow  coastal  waters  are  deep  enough  to  allow  the  operation 
of  a  fleet-class  submarine,  the  use  of  geophones  on  their  ice  covers  should 
be  considered.  In  thi.s  section  we  present  the  calculated  signal  intensities 
for  guided  waves  propagating  within  a  generalized  shallow  water  environment. 
This  we  have  taken  to  be  a  homogeneous  water  layer  100  meters  thick  underlain 
by  a  homogeneous  half-space.  The  characteristics  of  this  model  environment 
are  shown  in  Table  2.  Our  presentation  of  the  numerical  results  is  the  same 
as  employed  in  the  discussion  of  the  Arctic  Ocean  environment. 


24 


CONFIDENTIAL 


CONFIDENTIAL 


Table  2 

The  Approximation  for  the  Shallow  Water  Environment 


Layer 

Velocity 

Density 

Thickness 

(m/sec) 

(g/ cm^) 

(m) 

1436 

1,026 

2 

00 

For  the  environmental  parameters  listed  in  Table  2  the  phase  velocity 
dispersion  for  various  modes  was  calculated  and  is  shown  in  Fig.  12.  Based 
on  this  information  the  relative  signal  strengths  in  each  of  the  first  10 
modes  was  calculated  for  source  depths  of  25  and  50  meters,  Figs.  13  and  14, 
respectively.  These  strengths  were  then  combined  with  the  range  information 
in  a  manner  described  earlier  to  yield  the  variation  of  particle  velocity 
amplitude  with  range  at  frequencies  of  20,  50,  and  100  Hz.  We  again  assume 
unit  power  output  at  all  frequencies. 


FIG.  12  CALCULATED  PHASE  VELOCITY  DISPERSION  OF  GUIDED  WAVES 
IN  SHALLOW  ARCTIC  WATER 
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FIG.  13  EXCITATION  STRENGTH  OF  FIRST  10  NORMAL  MODES 

IN  SHALLOW  ARCTIC  WATER,  SOURCE  DEPTH  25  METERS 


FIG.  14  EXCITATION  STRENGTH  OF  FIRST  10  NORMAL  MODES 

IN  SHALLOW  ARCTIC  WATER,  SOURCE  DEPTH  50  METERS 
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For  a  20'Hz  source  the  amplitudes  are  plotted  against  range  in  Fig.  15. 

The  uppermost  curve  (dashed)  corresponds  to  the  signal  amplitudes  from  a 
source  at  a  depth  of  25  meters  when  the  roughness  of  the  ice  bottom  is  ig¬ 
nored.  However,  in  a  shallow  water  environment  the  number  of  reflections 
per  waveguide  distance  is  larger  than  in  a  deep  water  environment,  and  the 
scattering  of  incident  waves  by  the  ice  irregularities  causes  a  significant 
attenuation,  even  for  frequencies  as  low  as  20  Hz.  This  may  be  seen  in 
Fig.  15  by  comparing  the  dashed  curve  to  the  solid  curve  for  a  source  depth 
of  25  meters. 

We  do  not  know  if  we  are  justified  in  using  for  this  environment  the 
attenuation  factor  based  on  a  3-meter  standard  deviation  in  ice  roughnes  !. 

We  rather  suspect  that  this  value,  and  thus  the  resulting  attenuation,  is 
too  large  in  relation  to  the  smaller  ice  thicknesses  found  in  some  ice- 
covered  coastal  regions.  However,  since  we  do  not  include  the  effects  of 
other  causes  for  attenuation  in  shallow  water,  such  as  shear  wave  conversion 
and  bottom  irregularities,  the  correction  factor  probably  does  not  yield 
unrealistic  values  of  signal  amplitude. 

The  ice-bottom  roughness  results  in  an  unequal  attenuation  of  the 
modes  of  propagation.  This  is  indicated  in  Fig.  15  by  the  phase  interaction 
between  the  first  two  modes,  which  is  negligible  beyond  a  range  of  8  kilometers 
after  we  include  the  roughness  attenuation.  The  overall  rate  of  attenuation 
appears  to  be  slightly  greater  for  a  25-meter  source  depth  than  a  50-meter 
source  depth;  a  feature  thought  to  be  caused  by  the  disparity  in  mode  excita¬ 
tion  strrngths  with  changing  source  depth. 

To  our  knowledge  there  has  been  no  experimental  work  to  determine  the 
range  of  noise  amplitude  on  ice  in  a  shallow  water  environment;  in  Fig.  15 
we  have  included  for  illustration  only  the  maximum  and  minimum  noise  level 
of  deep  sea  Arctic  ice  at  20  Hz.  If  these  noise  levels  apply  equally  to  both 
acoustic  environments,  there  will  be  a  more  favorable  signal-to-noise  ratio 
in  the  shallow-water  case  than  at  a  comparable  frequency  in  deep  water.  This 
'uoulri  bo  the  case  if  more  of  the  noise  vibrations  at  20  Hz  are  generated 
locally  in  the  ice  and  propagate  principally  through  the  ice  layer.  On  the 
other  hand,  the  assumption  that  the  noise  is  propagated  through  the  water 
from  distant  sources  as  normal  mode  vibrations  would  imply  that  the  signal- 
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to-noise  ratio  in  the  shallow-water  case  could  be  as  poor  as  in  the  deep¬ 
water  case.  Tliat  is,  the  shallow  water  environment,  conducive  to  the  strong 
excitation  of  normal  mode  vibrations,  does  not  distinguish  between  the  prop¬ 
agation  of  signal  and  noise  disturbances.  However,  it  is  more  likely  that 
ice  vibration  noise  is  local  and  that  a  more  favorable  signal-to-noise  ratio 
occurs  in  the  shallow  water  case.  Hence,  despite  the  pronounced  attenuation 
caused  by  ice  roughness,  the  signal  amplitude  from  a  one-watt  source  could 
be  equal  to  or  greater  than  the  likely  noise  amplitude*  out  to  ranges  of 
14  kilometers. 

Figures  16  and  17  show  the  calculated  signal  intensities  at  50  and  100  Hz 
for  sources  25  and  50  meters  deep  after  the  correction  for  ice-bottom  scattering 
has  been  made.  The  scattering  causes  a  rapid  attenuation  of  signals  at  these 
higher  frequencies.  As  a  result  the  anticipated  signal-to-noise  ratios  would 
be  much  less  favorable  than  those  in  the  20-Hz  case.  A  comparison  between 
Figs.  11  and  17  reveals  that  the  average  signal  level  at  100  Hz  is  no  better 
in  shallow  water  than  in  deep  water  for  sources  of  equal  strength. 

Comparison  of  Hydrophone  and  Geophone  Signal-to-Noise  Ratios 

Comparisons  have  been  made  between  signal-to-noise  ratios  for  geophones 
and  hydrophones  equidistant  from  distant  underwater  explosive  sources  (Greene, 
Elbert,  and  Fitzpatrick;  1966).  Relative  signal-to-noise  ratios  are  shown 
in  Fig.  18,  taken  from  Greene  et  £l.  The  general  conclusion  from  this  work 
is  that  the  hydrophone  is  a  much  better  sensor,  regardless  of  frequency  and 
source  range.  This  is  perplexing  at  first  glance.  We  would  not  expect  the 
significant  difference  in  the  calculated  ratios  to  occur  because  both  pressure 
and  particle  velocity  are  analytically  related  to  the  same  acoustic  potential 
functions  that  describe  signal  and  noise.  Moreover,  the  relatively  poorer 
ratios  from  the  geophone  data  cannot  be  explained  away  by  signal  attenuation 
in  the  ice  at  the  geophone  site.  The  data  given  in  the  figure  do  not  suppox’t 
this  explanation. 


*We  assume  here  a  likely  noise  amplitude  of  65  mp,/sec,  a  value  explained  in 
the  previous  section. 
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FIG.  16  VERTICAL  PARTICLE  VELOCITY  AS  A  FUNCTION  OF  RANGE, 
SHALLOW  ARCTIC  WATER,  50-Hz  OUTPUT 
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FIG.  18  COMPARISON  OF  HYDROPHONE  TO  GEOPHONE  SIGNAL-TO-NOISE  RATIO 
FOR  SIGNALS  RECEIVED  FROM  DISTANT  EXPLOSIONS 
(After  Greene,  Elbert,  Fitzpatrick,  1966) 


It  is  our  belief  that  the  ice  is  a  relatively  noisier  environment  than 
IS  the  water  at  a  depth  of  one  hundred  feet  or  so  beneath  the  ice.  This 
state  could  occur  if  seismic  disturbances  originating  in  the  ice  are  effec¬ 
tively  trapped  within  the  ice  layer,  decaying  exponentially  in  amplitude  into 
the  water.  Various  types  of  trapped  waves  are  possible  in  a  floating  ice 
sheet  (Press  and  Ewing,  1951);  for  example,  the  ice  admits  trapped  longi¬ 
tudinal  waves.  This  hypothesis  does  not  explain  the  better  relative  geophone 
signal-to-noise  from  the  nearest  explosion  reported  (49  n.m.).  Perhaps  this 
was  due  to  quiet  ice  conditions  at  the  time  of  the  experiment. 

Over  the  20  to  100-Hz  bandwidth,  specifically  studied  in  this  report. 

Fig.  18  shows  tha'C  roughly  a  20  dB  better  signal-to-noise  ratio  is  obtained 
with  a  hydrophone.  The  data  further  indicate  that  above  100  Hz  the  hydro¬ 
phone  gives  a  40  dB  signal-to-ncise  improvement  over  the  geophone.  The 
marked  difference  above  100  Hz  implies  that  noise  sources  in  the  ice  must  be 
extremely  local  (the  ice  strongly  attenuates  high-frequency  vibrations). 
However,  as  we  do  not  know  the  frequency  response  of  the  geophone  nor  the  sys¬ 
tem  noise,  we  cannot  be  sure  that  the  40  dB  signal-to-noise  difference  is  real 
at  the  higher  frequencies  where  geophone  sensitivity  and  the  signal  strength 
fall  off. 
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Recognition  and  Localization  of  tie  Source 

Recognizing  and  localizing  a  nearby  underwater  noise  source  in  deep 
water  with  geophones  is  severely  restricted  by  the  natural  ice  vibrations, 
a  fact  brought  out  in  the  previous  sections.  Here  we  will  not  attempt  to 
deal  with  recognition  and  localization  in  terms  of  acoustical  outputs 
measured  from  submarines  of  various  classes  and  traveling  at  different 
speeds.  Our  purpose  is  only  to  bring  up  a  number  of  general  points  that 
have  a  bearing  on  the  usefulness  of  geophones  (and  hydrophones  as  well). 

The  usual  technique  for  improving  signal-to-noise  ratio  is  to  deploy 
numerous  sensors  in  arrays.  Earthquake  seismologists  are  well  accustomed 
to  the  idea  of  giant  arrays.  The  large-aperture  seismic  array  (LASA)  in 
Montana  consists  of  over  500  seismometers  in  21  clusters  over  a  circular 
area  of  200-km  diameter.  Recently,  it  was  suggested  (Press  and  Brace, 

1966)  that  about  1000  seismometers  be  arrayed  in  selected  locations  of 
California  and  Nevada  as  part  of  an  experiment  to  detect  and  locate  earth¬ 
quakes  and  micro-earthquakes.  When  the  signal  sought  is  coherent  over  a 
large  area  and  the  noise  is  not,  array  processing  of  seismometer  outputs 
has  led  to  an  increase  in  the  signal-to-noise  ratio  by  a  factor  slightly 
greater  than  the  square  root  of  the  number  of  sensors. 

However,  the  two  acoustic  environments  under  discussion  are  dispersive, 
and  this  results  in  signals  that  are  not,  in  general,  spatially  coherent. 

Tlie  unhappy  coincidence  of  unfavorable  signal-to-noise  ratio  and  dis¬ 
persive  signals  does  not,  in  theory,  eliminate  the  possibility  of  localizing 
the  source.  However,  a  large  number  of  on-ice  transducers  must  be  used  for 
this  purpose. 

In  a  homogeneous  medium  the  simplest  means  for  localizing  a  source  is 
a  steerable  array.  The  signals  sensed  by  several  colinear  transducers  are 
brought  into  phase  by  the  introduction  of  various  time  delays.  The  necessary 
delays  are  related  to  the  source  azimuth. 

In  an  inhomogeneous  medium  where  guided  wave  propagation  occurs,  each 
frequency  component  will  travel  with  a  multiplicity  of  phase  and  group 
velocities,  all  of  which  will  be  different.  In  this  situation  time  delay 
steering  will  not  yield  the  true  source  azimuth.  For  example.  Clay  (1961) 
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studied  array  steering  for  the  localization  of  a  continuous  signal  source 
in  a  shallow-water  waveguide,  and  showed  that  interference  effects  between 
modes  led  to  a  multiplicity  of  source  azimuths.  In  this  case  three  modes 
were  dominant  and  three  apparent  source  azimuths  were  indicated  by  the 
processed  hydrophone  output  signals:  one  apparent  azimuth  was  real  but 
wrong  and  two  apparent  azimuths  were  imaginary.  The  source  was  colinear 
with  the  hydrophone  array  in  Clay's  study. 

Only  in  the  special  case  where  the  source  azimuth  is  +  n/2  radians 
with  respect  to  the  direction  of  the  linear  array  may  the  source  azimuth 
be  correctly  determined.  In  that  case  the  traveling  wave  fronts  are  parallel 
to  the  line  of  detectors  and  the  signals  sensed  a.'e  coherent.  Tliis  .suggests 
that  if  we  were  to  deploy  numerous  linear  arrays,  each  like  a  spoke  radiating 
from  the  hub  of  a  wheel,  the  linear  element  whose  summed  and  squared  signal 
is  a  maximum  is  most  nearly  parallel  to  the  signal  wave  fronts. 

The  number  of  sensors  necessary  in  each  linear  element  is  a  subject  of 
concern.  The  more  sensors  used  the  greater  the  signal  enhancement  and  the 
less  chance  of  correlating  noise  and  coming  up  with  an  erroneous  source  azi¬ 
muth.  On  the  other  hand,  a  large  number  of  sensors  adds  to  the  operational 
problems.  For  localizing  a  source  emitting  one  watt  at  7  Hz  a  6  to  10  dB 
signal-to-noise  enhancement  would  be  desirable,  and  this  would  require  at 
least  four  on-ice  sensors  in  each  linear  element  of  the  array, 

The  only  benefit  trom  a  single  linear  array  is  indication  of  whether 
a  quasi-stationary  acoustic  source  is  present  in  the  vicinity.  The  square 
of  the  summed  and  delayed  sensor  outputs  is  given  by  two  terms  analytically, 
one  of  which  is  the  covariance  between  sensor  outputs.  As  Clay  (1961)  showed, 
in  a  dispersive  medium  the  time  delays  between  sensor  outputs  that  yield  a 
maximum  covariance  cannot  be  used  to  estimate  source  azimuth.  However,  if 
these  delay  increments  change  slowly  and  systematically  over  a  fairly  long 
period  of  time,  the  likelihood  of  a  quasi-stationary  source  (submarine) 
would  be  suggested. 

With  some  reservations  we  point  out  that  source  localization  in  a  noisy, 
dispersive  medium  is  possible  if  we  were  to  ignore  phase  relations  and  only 
the  signal  intensity  outputs  from  many  geophones  arrayed  over  a  wide  area 
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around  a  suspected  source.  It  was  illustrated  graphically  in  the  preceding 
sections  that  the  phase  velocity  dispersion  of  acoustic  signals  gives  rise 
to  a  signal  intensity  that  falls  off  with  range  in  an  oscillatory  manner. 

Thus,  the  signal  intensities  sensed  by  just  a  few  widely  spaced  geophones 
near  a  submarine  are  not  accurate  indicators  of  the  source  position;  a 
larger  signal  intensity  might  be  sensed  at  a  site  farther  from  tlie  source 
than  at  one  closer.  Moreover,  the  noise  intensity  and  the  non-uniformity 
of  noise  intensity  would  add  to  the  uncertainty.  However,  if  we  were  to 
"trap"  the  source  v/ithin  a  closely  spaced  network  of  many  sensors  we  could 
use  a  comparison  of  individual  sensor  outputs  to  identify  the  source  position. 
This  is  a  sledge  hammer  approach. 

The  approach  is  less  demanding  on  the  number  of  geophones  where  the 
source  position  is  constrained  by  natural  geographic  barriers — such  as  the 
ice-covered  waters  of  the  Canadian  Archipelago.  There,  acces.s  routes  for  a 
submarine  are  limited  in  places  to  deep  water  channels,  and  the  number  of 
geophones  becomes  the  square  root  of  the  number  required  to  yield  the  same 
coverage  in  an  unbounded  environment. 

The  Use  of  On-Ice  Sensors  Other  than  Geophones 

Although  the  use  of  geophones  has  been  emphasized  in  this  report,  we 
recognize  that  a  wide  variety  of  transducer  types  exist  for  the  detection 
of  small  ice  vibrations,  e.g. ,  piezoelectric  and  magnetostrictic  sensors. 

These  transducers  are  normally  used  to  measure  high-frequency  vibration 
(above  10  kHz)  whereas  we  have  focused  on  the  low-frequency  sensitive  geo¬ 
phone.  However,  in  view  of  the  unfavorable  ratio  predicted  between  the 
signal  from  ship-generated  noi.se  nnd  natural  ice  vibration  noise  at  low 
frequencies,  we  are  prompted  to  disc.iss  briefly  the  potential  use  of  the 
high-frequency  sensitive  transducers. 

A  high-frequency  sonar  signal  might  be  transmitted  from  the  submarine 
for  detecting  dangerous  under-ice  protuberances.  Because  these  sonar  signals 
are  above  1  kHz  the  total  rate  of  signal  attenuation  due  to  spreading,  absorp¬ 
tion,  refraction  into  ice,  and  scattering  from  the  rough  ice  under-surface 
will  be  much  greater  than  for  the  low  frequencies  previously  discussed. 
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However,  the  unfavorable  attenuation  is  offset  by  two  important  factors: 
the  natural  vibration  amplitude  '  crease  with  increasing  frequency,  and 
the  acoustic  power  output  of  the  sonar  transducer  is  much  larger  than  the 
power  in  spectral  lines  from  submarine  hydrodynamic,  flow,  and  machinery 
noise  sources. 

To  our  knowledge  natural  ice  vibration  studies  have  not  yet  been  made 
above  1  kHz.  However,  if  the  Knudsen  noise  curves  indicate  the  variation  of 
ice  noise  with  frequency,  then  the  noise  amplitudes  could  be  20  to  40  dB 
below  the  noise  level  for  the  10-  to  lOQ-Hz  bandwidth. 

The  potential  use  of  high-frequency  on-ice  transducers  is  an  intriguing 
subject  but  the  calculation  of  high  frequency  signal  vibrations  is  not  tractable 
to  a  simple  analytical  study  such  as  given  in  this  report:  the  short  wavelengths 
of  acoustic  energy  above  1  kHz  pose  greater  complexity  to  the  estimation  of 
vibration  amplitudes  than  at  lower  frequencies. 

Conclusions  and  Recommendations 

The  studies  reported  here  indicate  that  below  100  Hz  ambient  ice  noise 
in  the  central  Arctic  Ocean  results  in  a  generally  unfavorable  signal-to-noise 
state  for  the  on-ice  detection  of  a  nearby  submarine.  In  the  1  to  10-Hz  range 
the  noise  amplitudes  are  reported  to  be  a  minimum,  and  it  is  suggested  that  a 
submarine  emitting  a  significant  amount  of  power  (one-watt)  at  7  or  8  Hz  would 
produce,  within  a  20  or  25  km  radius  from  the  source,  a  vibrational  signal  in 
the  ice  roughly  equal  to  the  amplitude  of  the  maximum  noise  vibration.  From 
10  to  100  Hz  the  natural  ice  vibrations  are  uniformly  large  and  the  likely 
noise  amplitudes  are  sufficient  to  overwhelm  the  signal  from  a  nearby  source 
emitting  one-watt  at  discrete  frequencies. 

These  conclusions  are  based  on  numerically  derived  signal  amplitudes  for 
a  point  source  in  an  acoustical  environment  approximating  the  central  Arctic 
Ocean,  and  on  natural  ice  vibration  noise  measurements  reported  for  the  central 
Arctic  Ocean. 

The  noise  data  used  are  from  observations  made  during  the  late  winter 
months.  We  would  expect  the  natural  vibrations  to  be  ower  during  summer 
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months  when  the  ice  becomes  "fluidized"  and  a  poorer  generator  of  seismic- 
acoustic  noise,  and  thus  the  noise  data  used  represent  a  worst  noise  state 
for  the  central  Arctic  Ocean.  This  state  is  most  pertinent  to  the  study 
for  it  corresponds  to  the  season  when  the  ice  is  thickest  and  leads  fewest, 
consequently  when  hydrophone  use  is  most  severely  restricted. 

A  literature  search  showed  that  noise  vibrations  have  not  been  syste¬ 
matically  measured  in  Arctic  ice  that  covers  shallow-water  regions.  Tins 
is  unfortunate  for  the  noise  spectra  may  not  be  identical  in  both  deep- 
and  shallow-water  regions.  Nevertheless,  signal  amplitudes  were  also  cal¬ 
culated  for  point  sources  in  an  idealized  shallow-water  acoustical  environment 
with  a  water  depth  of  100  meters.  At  frequencies  below  about  50  Hz  the  cal¬ 
culated  signal  amplitudes  at  short  ranges  (<  15  km)  are  as  much  as  20  dB 
larger  in  shallow  water  than  in  deep  water.  This  difference  is  caused  by 
the  stronger  excitation  of  normal  mode  vibrations  in  the  shallow  water. 
However,  above  50  Hz  the  ice-roughness  attenuation  of  acoustic  waves  in 
shallow-water  regions  offsets  the  stronger  excitation  of  the  waves. 

If  the  ice  vibration  noise  above  5  Hz  is  caused  principally  by  local 
disturbances  occurring  in  and  propagating  through  the  ice,  than  the  noise 
levels  could  be  roughly  the  same  in  both  environments.  With  this  assumption 
the  signal-to-noise  ratio  from  5  to  50  Hz  might  be  20  dB  better  at  short 
ranges  in  shallow  water  than  in  the  central  Arctic  Ocean.  Thus,  the  on-ice 
transducer  which  seems  to  have  at  best  marginal  value  in  deep-water  detection 
might  be  extremely  effective  for  shallow-water  detection  of  submarines. 

Whether  this  is  true  can  only  be  determined  from  future  ice  vibration  studies 
in  shallow-water  regions  of  the  Arctic. 

From  tests  in  the  central  Arctic  Ocean  it  has  been  shown  that  the  hydro¬ 
phone  is  much  more  effective  in  detecting  underwater  sound  sources  than  the 
geophone.  Comparisons  of  signals  observed  by  menus  of  the  two  transducers 
show  a  10  to  40  dB  signal-to-noise  improvement  with  the  hydrophone.  This  is 
explained  by  the  noise  waveguide  behavior  of  the  ice  layer:  the  ice  traps 
much  of  the  vibrational  energy  generated  within  it.  Consequently,  a  hydro¬ 
phone  must  be  preferred  over  an  on-ice  transducer  whenever  an  alternative 
exists. 
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Wlien  there  is  no  choice  but  to  use  an  on-ice  transducer  the  transducer 
used  should  have  maximum  sensitivity  over  the  bandwidth  of  maximum  signal- 
to-noise  ratio.  Ambient  noise  is  a  minimum  from  2  to  10  Hz,  and  submarine 
generated  acoustic  signals  might  be  relatively  large  in  this  frequency  range. 
Therefore,  a  geophone  with  a  resonant  frequency  from  2.0  to  7.5  Hz  is  sug¬ 
gested  for  on-ice  sensing  in  deep  Arctic  waters. 

Despite  the  generally  unfavorable  signal-to-noise  ratio  and  dispersion 
of  the  signals,  geophones  can  theoretically  be  used  in  two-dimensional  arrays 
for  estimating  source  azimuth.  However,  since  one  hundred  or  more  transducers 
are  necessary  the  operational  problems  of  deploying,  monitoring  and  processing 
the  geophone  outputs  would  be  severe.  Essentially,  the  same  can  be  said  for 
hydrophone  operation.  Although  the  signal-to-no,lse  ratio  of  the  hydrophone 
output  should  be  better,  the  signal  sought  is  still  dispersive  and  two- 
dimensional  hydrophone  arrays  are  necessary  for  source  localization. 

Existing  data  concerned  with  natural  ice  noise  is  grossly  inadequate. 
Experimental  work  as  follows  is  therefore  recommended. 

Ice  noise  should  be  measured  on  ice  over  water  less  than  100  meters 

deep. 

The  frequency  range  should  be  from  1  Hz  to  the  frequencies  of  submarine 
sonar.  However,  the  emphasis  should  be  on  the  frequencies  between  1  Hz  and 
50  Hz. 

These  measurements,  particularly  the  low  frequency  ones,  should  be  made 
with  three-component  geophones.  Sound  waves  from  sources  within  the  ice  layer 
are  likely  to  have  important  horizontal  components  as  well  as  vertical  com¬ 
ponents,  whereas  waves  originating  in  the  water  are  likely  to  have  primarily 
vertical  competent s  in  the  ice. 

Measurements  of  ice  noise  should  be  carried  out  to  the  extent  that  some 
statistical  data  are  developed  regarding  the  variation  of  ice  noise  from  quiet 
days  to  stormy  days. 
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APPENDIX  A 

Determination  of  the  Eigen  Wavenumbers 

It  was  shown  by  Tolstoy  (1956,  1958)  that  the  characteristic  equation 
for  the  eigen  wavenumbers  of  totally  reflected  guided  waves  could  be  ex¬ 
pressed  in  terms  of  phase  angles  )(t  and  for  the  up-going  and  down-going 
waves.  If  the  acoustical  environment  is  approximated  by  n  homogeneous  layers, 
the  characteristic  equation  is 

■  e~‘X‘-  =  1 

or 


where  m  is  an  integer,  the  mode  number;  k  is  the  vertical  component  of 

mz 

the  eigen  wavenumber,  the  parameter  that  also  satisfies  the  phase  integral; 
and  hi  is  the  thickness  of  the  low  velocity  layer.  If  the  upper  boundary 
is  a  free  surface,  then  =  0  and  the  characteristic  equation  reduces  to 

=  mrr  (A.l) 

To  solve  for  the  eigen  wavenumbers  we  must  expand  terms  of  the  entire 

acoustical  system.  In  the  Arctic  Ocean  the  sonic  velocity,  a,  increases 

monotonically  with  depth.  Partitioning  the  environment  into  n  layers,  and 

assuming  no  partial  reflections,  total  reflection  will  occur  at  the  i, 

(i  +  1)  is  interface  where  the  phase  velocity  of  the  guided  wave,  c, 

a  <c<a  (o',  and  a.  ,  are  the  acoustic  velocities  in  the  (i  +  1)  layers, 
i  i+1  1  1+1 

respectively).  At  this  interface  the  resulting  phase  change  is 
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^1 ,  1  +  1 


where 


and  p  ,  p.  are  the  densities  of  the  i  and  (i  +  1)  layers,  respectively 
jL  i  4'X 

(Brekhovskikh,  1960).  In  the  simplest  of  all  cases,  a  single  low  velocity 
layer  over  a  half  space  (i  =  1),  we  find  that  for  a  given 


we  may  solve  explicitly  for  the  k 


mz 


However,  when  we  consider  a  less 


trivial  case  in  which  the  down-going  waves  penetrate  more  than  one  layer, 

the  determination  of  k  becomes  less  direct.  Each  layer  above  the  reflec 

mz 

ting  interface  will  affect  the  total  phase  change  experienced  by  the  wave. 
In  this  case  >3®  found  through  the  recursive  relationships  given: 


=  tan 


- 1 


Po 


X23  " 


Pz 

P23-—  tan  (Tj/ij  +  Xii) 
Po 


-  1  ,  n 


tan 


Pn 

Pn-l.n  ^  tan  (7„h„  +  ,  „  + ,  ) 

^n-1 


Y  4. ,  =  tan 

''■n  ,  n  +  1 


- 1 


P  n  ,  n  +  I 


^,:  +  ! 


Pn 


(A. 2) 
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where  /  is  the  magnitude  of  the  eigen  wavenumber  in  the  i  -  th  la\er, 


To  solve  for  Tolstoy  (1955)  gives  an  iterative  method  easily  adapted 

to  high-speed  machine  evaluation. 
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APPENDIX  B 

The  Effect  of  Ice  Roughness  on  the  Guided  Waves 

If  we  consider  a  harmonic  wave  with  amplitude  A  and  wavenumber  k  * 

m  mz 

incident  at  a  rough  or  irregular  sea  water-sea  ice  boundary,  then  the 
reflection  will  not  be  entirely  specular.  A  portion  of  the  incident  energy 
will  be  reflected  in  directions  whose  corresponding  wavenumbers  no  longer 
satisfy  the  condition  for  unattenuated  guided  waves.  This  energy  is  not 
trapped  in  the  guide  and  attenuates  rapidly.  Assuming  the  ice  thickness 
to  be  small  compared  to  the  wavelength  of  the  incident  wave  and  the  surface 
to  have  a  Gaussian  distribution  function,  then  total  reflection  of  the  in¬ 
cident  wave  occurs;  but  the  amplitude  specularly  reflected  portion  is  (Clay, 
1964) 


A' 


A.e 


where  the  random  surface  ^(x)  is  characterized  by  the  distribution  function 

D(i)  =  - ;  exp  (-^-/2cr2) 

cr(277)^^ 


Because  k  is  dependent  on  the  mode,  different  modes  will  suffer  different 
mz 

rates  of  attenuation.  In  general  k  increases  with  mode,  and  hence  the 

mz 

lowest  modes  should  contribute  most  to  the  signal  at  long  ranges. 

In  a  waveguide  distance  R  the  wave  will  suffer  n  reflections.  If  the 
thickness  of  the  guide  is  h,  the  number  of  reflections  is 


n  -  R/2h  tan  0 


*As  before  we  define  k^^  as  the  vertical  component  of  the  eigen  wavenumber. 
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where  0  is  the  angle  of  incidence  for  the  m-th  mode.  Thus  if  A  is  the 
m  ^ 

unattenuated  amplitude,  the  attenuated  amplitude,  (R)  after  waveguide 
distance  travel  R  becomes 


A'iB) 

A 


-2n  (  t 

e 


a  i 


■) 

)  - 


or 


~y . 


ft 


where 


and 


k-  * 

n  t 


When  more  than  a  single  mode  Is  involved  the  ice  roughness  will  also 
modify  the  signal  intensity  at  distance  R  by  causing  an  incoherence  of  the 
modes  relative  to  each  other.  In  Eq.  7  it  was  shown  that  the  signal  in¬ 
tensity  is  derived  in  part  from  a  term  expressing  the  phase  interaction  of 
the  modes.  This  term  has  a  dependence  on  range  R  like 

cos  {k  -  k 

nr  nr 

However,  in  a  waveguide  of  changing  thickness  neither  k^^  nor  k^^  are  constant, 

these  eigenvalues  change  slightly  with  perturbations  in  the  height  of  the  ice 

bottom  irregularities.  Clay  (1964)  states  that  cos(k  -  k  )R  must  be  re- 

mr  nr 

placed  by  a  term  derived  from  averaging  cos(k  -  k  )R  over  the  Gaussian 

mr  mz 

distribution,  D(S)>  Thus,  the  phase  interaction  term  is  shown  to  vary  with 
range  like 

(exp  {-(k  -  k  cos  (k  -  k  )R 

*  fflrnr  m  r  n  r 
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